Spin Injection in Quantum Wells with Spatially Dependent Rashba Interaction 
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We consider Rashba spin-orbit effects on spin transport driven by an electric field in semiconductor 
quantum wells. We derive spin diffusion equations that are valid when the mean free path and the 
Rashba spin-orbit interaction vary on length scales larger than the mean free path in the weak spin- 
orbit coupling limit. From these general diffusion equations, we derive boundary conditions between 
regions of different spin-orbit couplings. We show that spin injection is feasible when the electric 
field is perpendicular to the boundary between two regions. When the electric field is parallel to 
the boundary, spin injection only occurs when the mean free path changes within the boundary, in 
agreement with the recent work by Tserkovnyak et al. cond-mat/0610190]. 
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I. INTRODUCTION 

Spintronics envisions electronic devices with spin in- 
jection, detection, and manipulation. The spin-orbit 
interaction couples the electron spin to its orbital mo- 
tion. It was early realized that it can be utilized to 
manipulate the electron spin in e.g. the Datta-Das spin 
transistor, 1 - More recently, the possibility of spin injection 
via the spin-orbit interaction in non-magnetic systems 
was suggested^ 3 - The spin-orbit interaction gives rise to 
the spin Hall effect where a longitudinal electric field in- 
duces a transverse spin current. In semiconductors, the 
spin Hall effect has generated a large interest and is pre- 
dicted to occur in a wide variety of electron and hole 
doped 2D and 3D systems^ A resulting accumulation of 
a normal to 2DEG spin polarization near sample bound- 
aries has recently been measured^ Also, injection of the 
spin Hall polarization into a region of the 2DEG where 
the driving electric field is absent has been observed in 
Ref. S 

The spin Hall effect is conventionally separated into an 
intrinsic and an extrinsic effect. The intrinsic spin Hall 
effect is caused by the spin-orbit splitting of energy bands 
in zinc-blende semiconductors. In contrast, the extrinsic 
spin Hall effect arises due to spin-dependent scattering off 
impurities. This paper considers the intrinsic spin Hall 
effect in the diffusive transport regime for electrons con- 
fined in a quantum well subject to the Rashba spin-orbit 
interaction. We will focus on spin diffusion in systems 
where the Rashba coupling constant a varies in space. In 
particular, we address a problem of injection of the spin 
polarization from regions with a large spin-orbit coupling 
a to regions where this constant is small or zero. 

Recent research has established how spins diffuse, pre- 
cess and couple to the charge transport in electron doped 
III-V semiconductor quantum wells. The charge and 
spin flows are described by spin-charge coupled diffusion 
equations;^ However, the diffusion equations, valid in 
the bulk of the systems, must be supplemented by ap- 
propriate boundary conditions. Calculations to this end 



were carried out in Refs. I^ ITfllTTIfl2lll3lll4l and have given 
rise to discussion in recent literature. The derivation of 
these boundary conditions is more subtle than it might 
first seem since the calculations must be carried out to 
the second order in the spin-orbit interaction strength 
and include the effects of the electric field. 

The correct boundary condition when the electric field 
is parallel to the boundary and th e sp in-orbit coupling 
varies in space was derived in Ref. 1 1 71 - An argument in 
terms of conventional Hall physics that extends the re- 
sults to systems where also the mobility varies was also 
given: A spatially dependent spin-orbit coupling varying 
on length scales smaller than the spin-precession length 
corresponds to an effective magnetic field of opposite 
magnitude for two spin directions ! 17 ' 19 The effective mag- 
netic field induces spin dependent Hall voltage across the 
boundaryji 7 - This additional voltage, in its turn, gives rise 
to the spin density variation in the transition layer, re- 
sulting in a finite spin-density jump, when the transition 
layer thickness is reduced to zero. 

This paper aims to address spin diffusion and the 
boundary conditions in more detail and generality for 
the Rashba spin-orbit interaction model. Our aim is 
to clarify some conceptual issues and not to fit exper- 
imental data because typical experimental systems are 
rather complex with competing effects and material pa- 
rameters that are not well known. The Rashba spin-orbit 
interaction in semiconductor quantum wells depends on 
quantum well confinement, which can be manipulated by 
gates, or doping profile^ but usually only at rather long 
scales. The length scales for variations in system prop- 
erties such as the spin-orbit interaction and the mean 
free path are typically much larger than the mean free 
path. Therefore, we will derive general diffusion equa- 
tions that are valid when the system parameters vary on 
length scales longer than the mean free path. These diffu- 
sion equations can then be solved for smooth boundaries 
between regions of different spin-orbit coupling strengths 
and mean free path. From these solutions, we will find 
boundary conditions that are not limited to a specific 
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relative direction between the electric field and the inter- 
face. When the electric field is parallel to the interface 
our derived boundary conditions agree with the result by 
Tserkovnyak et al^ and the analog with conventional 
Hall physics is rigorously justified beyond the simplest 
scenario considered microscopically in Ref. [TtI We also 
present results for the boundary conditions when the elec- 
tric field is perpendicular to the interface. 

Our paper is organized in the following way. In the 
next section UH we summarize our main results for the 
resulting boundary conditions between regions with dif- 
ferent spin-orbit couplings and mean free path. We dis- 
cuss the implications of these boundary conditions on 
spin injection and show that for a parallel electric field 
spin injection is only possible when the mean free path 
together with a vary within the boundary layer, while 
for the perpendicular field spin injection takes place even 
when the mean free path is uniform in the system. Using 
the Keldysh formalism, we derive the general diffusion 
equation in section IIIII that we employ to compute the 
boundary conditions presented in section [TXJ Our conclu- 
sions are in section IIVI 



II. 



BOUNDARY CONDITIONS AND SPIN 
INJECTION 



In the absence of disorder, the electrons in a homoge- 
neous quantum well are described by the Hamiltonian 
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(1) 



where m* is the electron effective mass, k = (k x ,k y ) is 
the wave vector, cr is a vector of Pauli matrices, and hk 
is the effective spin-orbit field. Let us assume that the 
major contribution to the spin-orbit interaction (SOI) is 
given by the Rashba interaction: 



ak v 



K 



-ak x 



(2) 



where ki is the 2D electron wave vector. We also as- 
sume that the spin-orbit coupling constant a(r) is spa- 
tially modulated. In this case, the spin orbit term in Eq. 
((T|) has to be modified to preserve the Hermitian form of 
the Hamiltonian, namely, this term is rewritten as 1/2 of 
the anticommutator of the momentum operator and the 
spatially dependent a(r). 

Since the spatial variations of a are assumed to be 
provided by an appropriate quantum well modulation or 
impurity doping profile, any consistent model must in 
general also take into account spatial variations of the 
mean free scattering time r(r). The electron gas is under 
the action of the electric field E(r) and our goal is to find 
the spin density induced in the 2DEG by E. The electric 
field may vary in regions where the mobility changes in 
order to maintain a constant charge current density. We 
will derive the general diffusion equations corresponding 
to the Hamiltonian (Q]) in the next section (JTTTJ) . 



Let us first discuss some of our main findings of the 
general diffusion equations when applied to the problem 
of spin injection from a region with a large spin-orbit 
coupling to a region with a smaller spin-orbit coupling. 
As an example, we consider a left (x < —d/2) and right 
region (x > d/2) with different spin-orbit couplings and 
mean free paths connected via a smooth transition layer, 
as shown in Figure [T] We assume that the transition 
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FIG. 1: A 2DEG with a spatial dependent Rashba interaction 
and scattering time. The Rashba interaction a(x) and scat- 
tering time t(x) vary in the transition region between the left 
region where a — ai and t — n and the right region where 
a = olt and r = r r . 

layer thickness d between the regions is larger than the 
mean free path I, but smaller than the spin precession 
length A = ti 2 /m*a, where m* is the electron effective 
mass, I <C d <C A. This requires a sufficiently weak spin- 
orbit interaction or a low mobility sample. We consider 
the general case where the electric field is in the two 
dimensional plane, E = E x x + E v y. 

The diffusion equations for the spin-components S x , 
S v , and S z of the spin accumulation in the bulk of each 
region (x < —d/2 or x > d/2) valid when 1/r, h^ F <C Ep 
are 
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DV 2 S X + ADm*aVS z 
DV 2 S y 

DV 2 S Z - ADm*aVS x 



= r„ (s* - sf) , 
= r l{ (sy-s y b ), 
= r±s z , 



(3) 



where D — v 2 f t/2 is the diffusion constant in terms of 
the Fermi velocity vp and the elastic scattering time r 
and r ± = 2F|| = 8D (m*) 2 a 2 are the D'yakonov-Perel 
spin relaxation rates. It is well known that in systems 
with Rashba interaction, a uniform electric field induces 
a noncquilibrium bulk spin polarization parallel to the 



and S^. 



Within the 



2DEGj±£ denoted in Eqs. © as S% 
linear response theory, the spin polarization components 
are 



S x = -aeE y TN F , 
qy _ 



aeE x TN F , 



where A^p = m* / (f\ 2 2ii) is the density of states at the 
Fermi energy. The spin current is not a conserved quan- 
tity in spin-orbit coupled systems, but, with the con- 
ventional definition within the semiclassical theory. 
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reads in diffusive systems for spin flow, e.g., along the x 
direction^ 

f = D [V X S X + 2m*aS z ] , 

f = D [V X S Z - 2m* a (S x - S x )} . 

We generalize below the diffusion equations to sys- 
tems with spatially varying parameters a and r. 

By solving these generalized diffusion equations in the 
boundary transition layer (—d/2 < x < d/2), that is thin- 
ner than the spin-precession length A, we find in Sec. Ill 
the boundary conditions 

f d/2 x^dr 

(S - S x ) d j 2 - (S x - S b )_ d j 2 = -J dxS b ~-,(4) 

^d/2 ~ S-d/2 = °> ' ( 5 ) 
$d/2 ~ S-d/2 = . (6) 

The boundary condition for the S x spin density compo- 
nent can be understood in terms of a conventional Hall 
effect for two spin directions in the transition layer, see 
the discussion in Rcf. We outline in the next section 
lllll how the boundary condition for the S v component can 
be understood as a result of electron drift in an effective 
spin dependent potential. 

The boundary condition for the spatial derivative of 
the spin density is particularly simple. It can be written 
as continuity of spin current in its conventional definition: 

(3s) d/ 2 ~ CU-d/2 = (7) 

Naively, the continuity of the spin current across the 
boundary is not surprising because the boundary con- 
ditions have been obtained assuming that the boundary 
layer is thinner than the spin precession length. That 
means that effects violating spin current conservation 
only occur on much larger length scales. Consequently, 
there is spin current conservation through the boundary. 
However, the continuity of the spin current across the 
boundary cannot be rigorously established without do- 
ing an actual calculation. There could in general be a 
violation of the spin current conservation of the order 
a 2 , but we find explicitly in the next section UHl that also 
such corrections are absent. 

Let us now discuss the implications of the diffusion 
equations ^ and the boundary conditions (JH), (0, ©, 
and ([7]) for spin injection from a region with e.g. a larger 
spin-orbit coupling constant to a region with a smaller 
spin-orbit coupling constant. We consider two limits for 
the alignment between the electric field and the bound- 
ary: 1) The electric field is parallel to the boundary, 
E = E v y, E x = and 2) the electric field is perpendicu- 
lar to the boundary, E = E x x, E y — 0. Since we consider 
the linear response regime, spin injection for general ori- 
entations between the electric field and the boundary can 
be found as a linear combination of 1) and 2). 



The former case 1) of a parallel electric field E = E v y 
corresponds to the case studied by Tserkovnyak et al<^ 
For this geometry, the bulk nonequilibrium spin accumu- 
lation along the y direction vanishes, S b — 0, and the 
boundary condition ([5]) implies that 

S v = 

throughout the system. The only nonvanishing bulk com- 
ponent of the spin distribution is along x. Deep inside 
the bulk of the left (1) and right (r) regions: 

S x t = -aieEV n N F , 
S x r = -a r eE y T r N F . 

If the electron mobility is the same in each region, n = 
iV, and constant throughout the boundary, then the spin 
accumulation S x will change rapidly (i.e., on the scale of 
the spatial variations of a, d) from S x = —aieE v TiN F 
at x = -d/2 to S x r = -a r eEyT r N F at x = d/2. For 
the slowly varying S x in bulk, this looks as a disconti- 
nuity in the spin accumulation at the boundary. When 
the mobility differs in the two regions, the magnitude of 
discontinuity in the spin accumulation across the bound- 
ary is modified according to Eq. ([?]). This result differs 
from Ref. [9|, where a continuous spin density across the 
boundary was assumed. Solving the diffusion equations 
©, the solution in the left region is (x < —d/2) 7 

S x = S x bi + J2 Sl n exp %n (x + d/2) , 

71=1,2 

SZ = E Sl, n exp + d/2) 

n=l,2 

and in the right region (x > d/2) 

S x - Sl+ ]T Sl^expq^i-x + d/2) , 

71=1,2 

S z = Sl. n expq r , n (-x + d/2), 

71=1,2 

where q n is the solution with positive real part of q n = 
(1/A)y 2(— 1 ± iV7). Continuity of the spin current, 

Eq. (O, combined with the boundary conditions for S x 
and S z , Eqs. g]) and ©, imply 

Si r ,n = ^ra^S^ , 
^i u n = *n,n^S ) 

where 

f d/2 1 dr 

AS X = — / dxS x (x)- — 
J -d/2 t dx 

and >ff , , xf and >c z are lengthy algebraic dimension- 
less coefficients that depend on qi >n , q r _ n and the diffusion 
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coefficients Di and D r . From these boundary conditions, 
it should be clear that spin injection, which we define as a 
spin accumulation that differs from the bulk spin density 
Sfc, i.e., Si ^ 0, is only feasible when the mobility varies 
within the boundary layer, since otherwise AS X vanishes. 
If the mobility is constant throughout the system, there 
is no spin injection from the left region to the right region 
or vice versa. In general, AS X depends on the specifics of 
the variation of the mobility and the spin-orbit coupling 
in the boundary layer. 

Next, let us study the second scenario 2) where the 
electric field is perpendicular to the boundary, E = E x x, 
E y = 0. In this case, we will find that spin injection 
is feasible even when the mobility is uniform throughout 
the system. In the bulk of the left (right) region, we have 
S y = Si = aieEfnN F {S v = a r eE x T r N F ), while the 
other spin components vanish. Since the electric current 
is continous across the interface, the electric fields on 
both sides of the boundary must differ in order to adjust 
to the different mobilities so that E x ti — E x r r . The 
boundary conditions (HJ), and (O then imply that 

S x = and S z = 

throughout the system. In contrast to the previous sce- 
nario 1), where the in-plane spin component S x exhibits a 
discontinuity across the transition layer when the mobil- 
ity is uniform, S y is continous here. Solving the diffusion 
equations ©, we find in the left region (x < —d/2) 



sy = si + s%exp\i 1 ( x - 



■ d/2) , 



the mobility is constant throughout the system. For the 
simplest case, where both the mean free path and the 
spin-orbit coupling linearly change within the boundary 
layer: 



we find 



a(x) = 
t(x) = 

AS X = 



■ + ai x 

— + - d {a r - ai ) 

+ n x 

— + -d {Tr - Tl) ' 



Oil 



N F eE y (r r - ti) 



and 



K-SJO =-(«r- a t ) Nj 



eE x r r + eEfn 



(8) 



For this simplest system of a linear variation of the spin- 
orbit coupling and mobility and assuming that one uses 
the same electric fields in the two scenarios, then 2) is 
more effective than 1) provided 



ai 



Oil + 0>r 



> 



n 



n 



e.g., when the change in the spin-orbit coupling constant 
is larger than the change in mobility between the regions. 
Here, we are crudely characterizing the spin-injection ef- 
ficiency by the jump in spin density with respect to the 
bulk levels. 



and in the right region (x > d/2) 

S y = Si + Sf exp X; 1 (-x + d/2) . 



III. DERIVATION OF DIFFUSION EQUATION 



Continuity of the spin current and the boundary cond- 
tions for S y imply 



en 

Sf 



(si 
K 



si)] 
si)} 



where 



— 



DrX- 1 



DiXr 1 + DrXV 1 



DiX~ L + D r Xr 



In the limit of weak spin-flip relaxation in the left region, 
A; > X r Di/D r , we find xf = -1, x y = and x z = 0. 
In contrast to scenario 1), spin injection is feasible even 
when the mobility is constant throughout the system. 
This means that spin injection is feasible independent on 
the spatial variation of mobilities. 

Let us now compare the spin-injection efficiency in sce- 
nario 1) (electric field parallel to the interface) and 2) 
(electric field perpendicular to the interface). It is clear 
that geometry 2) is more effective than geometry 1) when 



We will in this section derive general diffusion equa- 
tions that are valid when the mean free path and spin- 
orbit coupling vary on length scales longer than the mean 
free path. It is well known that in a homogeneous sys- 
tem an electric field induces a spin polarization parallel 
to the 2DEGt^ In an inhomogeneous system, an inter- 
play of this effect and the spin Hall effect can give rise 
to a more complicated spin distribution. A convenient 
tool to study such a system is the spin diffusion equa- 
tion. We start from the Keldysh formalism and obtain 
the Boltzmann equation for the spin distribution func- 
tion. In terms of the Boltzmann function, the spin den- 
sity is defined as S = X)k9k, and the nonequilibrium 
charge density as 2^ k <7^ h . The linearized equation for 
the nonequilibrium part of the Boltzmann function can 
be written in the form (see e.g. Refs. Il6lfr7t ) 



vS7gl - a(z x V) l3 f 4 



9k 



VftL + eE 



dk 



2 (gk x hk) 



9k 



(9) 



where v=k/m*, i 

g ( k° } = -h k S(Ek - 



= x,y,z, S E = 5(E k - E F )S/N F , 
Ep) is the equilibrium Boltzmann 
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function for spin, and z is the unit vector perpendicu- 
lar to 2DEG. The first two terms on the left-hand side of 
Eq. ([9]) describe variations of the Boltzmann function due 
to particle motion with spin-dependent velocity, the 3rd 
term is the spin precession in the Rashba field, the 4th 
term is the spin-dependent acceleration in the spatially 
dependent Rashba field, and 5th term is associated with 
the driving electric field. The right hand side of Eq. ^ 
represents collisions in the self-consistent Born approxi- 
mation caused by weak and isotropic scalar disorder. 

In the leading approximation, the charge component 
of the Boltzmann function is determined by the drift of 
electrons in the external field: 

kE 

gt = er—5(E k - E F ) (10) 

m* 

which can be spatially dependent through r = r(r). Due 
to this dependence, the second term on the left hand side 
of Eq. © becomes important in the following calculation 
of the spin polarization. We note that this takes place 
only for the electric field parallel to the interface. For 
a perpendicular field, as we have mentioned above, the 
combination tE x is continous across the boundary. 

The next step is to derive from ([9]) the diffusion equa- 
tion for the spin polarization. There are two ways to 
do that, depending on how fast r(r) and a(r) vary in 
space. If the corresponding length scale is shorter than 
the electron mean free path I, the diffusion equation can 
be applied only within regions where variations of r(r) 
and a(r) are slow, while the regions of their fast modula- 
tion have to be described by the original kinetic equation 
([9|) . For example, one can consider a thin transition layer 
between two regions characterized by noncqual uniform 
values of r(r) and a(r). In this case, Eq. © can be used 
to derive boundary conditions for diffusion equations on 
both sides of the interface. If variations of the parame- 
ters take place on length scales larger than I, the general 
diffusion equation valid in the entire 2D system can be 
straightforwardly derived from Eq. ([9]). We will follow 
the second procedure assuming slow enough variations of 
r(r) and a(r), throughout the system. We will see in 

I 



where r' = dr/dx. 

Let us consider the limiting case when the mobility 
through the scattering time r is homogenous, r' = 0. It 
can then be seen that the solution of Eqs. (fT5|) and (fT7|) 
is S z — 0, and S x = S£ . Since S% depends on x through 



the following that the characteristic spatial scale of the 
spin density variations dictated by the diffusion equation 
is the spin precession length A = h 2 jra*a. Hence, this 
length must be much larger than L which means also that 
hk <C 1/t. Then, following Ref. [l6|, the diffusion equa- 
tion can be directly obtained from Eq. (J9j> . To this end, 
let us denote 

oi Q ch Bn^ 

= ^ + «(z x V)^ + |^-V/4 - eE^|- . (11) 

The function can now be found by applying the oper- 
ator 

A" 1 « t [1 - («V)t + («V)r(uV)r] . (12) 

which is the inverse to A = i>V + r _1 expanded to the 
second order with respect to tvV, where v is the particle 
velocity. In this way, we obtain 

fifk = A" 1 [Q - 2/i k (<7 k x «k)] , (13) 

where rat = h^/h^, while to order a 2 the vector product 
in (fT3"|) can be expressed as 

g k xn k = A- 1 (Qxn k ) + 2r 2 /i k Q 1 , (14) 

where Q± — Q — «k(«kQ)- In the h 2 . terms we disre- 
garded gradient corrections. Furthermore, the A -1 oper- 
ator in Eqs. (fl"4"|) and (fT3")) can be expanded according to 
Eq. (fT2")l . Below, we will keep the terms up to order V 2 
and a 3 in expansions in V and a, as well as up to a 2 V in 
cross products. Substituting Eq. (fT4")) and Eq. fTi"]) into 
Eq. (|13p and taking a sum over k, we arrive at a closed 
diffusion equation for the spin density. For simplicity, 
we consider the case when a and r depend only on the 
x coordinate, while the direction of the electric field is 
arbitrary. This is sufficient to discuss, e.g., a boundary 
between two regions with different Rashba couplings. Fi- 
nally, the system of diffusion equations takes the form 



-Y^-St) = (15) 

-r„(s»-s*) = o (i6) 

-S X )-T ± S Z - 0, (17) 



a(x), this spin density component simply follows the spa- 
tial dependence of the Rashba coupling constant, as it 
was discussed in the previous section. The behavior of 
S v is different. If, for example, a varies from some finite 
value to zero across the transition layer, Eq. (fT6|) shows 



t' 

V X DV X (S X - S x ) + V X DS%— + 2Dm*aV x S z + 2V x Dm*aS z 

T 

\7 x D\7 x S y 

V X DV X S Z - 2Dm*aS x - - 2Dm*aV x (S x - S£) - 2V x Dm*a{S x 

T 

I 
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that the electric field perpendicular to the interface, E x , 
drives injection of S v into a region of e.g. a vanishing 
Rashba spin-orbit interaction, corresponding to scenario 
2) analyzed in detail in section HT1 

Let us consider a thin transition layer where a(x) and 
t(x) vary between their constant values on the left and 
on the right from this layer. In this case, the spin polar- 
ization is expected to change fast within the layer, while 
it varies much more slower outside the layer, if its width 
d is much less than the spin-precession length A. The 
magnitude of the spin density variation across the layer 
can be found directly from the diffusion equations (| 1 5|) - 
(|17p . Indeed, in the leading approximation with respect 
to d/X, when x is within the layer, we retain in the dif- 
fusion equations only the leading terms proportional to 
V 2 , . After integration of the first of the equations from a 
point x\ just to the left from the interface to some point 
x, it becomes 

Va (S» - St) + sf- - [V X (S* - Sf) + S x b -] X=X1 . (18) 

T T 

The expression in square brackets can be set to 
because V x S%= Xl ~ (S x )/X, V X S% \ x = Xt = 0, and 
Sfr' /t\ x — Xi = 0. Hence, integration of Eq. (fT8|) from 
X\ to X2 yields 

(S x -St)\ x=X2 -(S x -SZ)\ x = Xx = - J dxS?~. (19) 

Similarly, it is easy to show that the z spin component 
slowly varies across the boundary, while the y spin com- 
ponent changes according to 

s y \ x=X2 -sy\ x=xl =o. (20) 

Equation (fT9f coincides with the result obtained in Ref. 
[YR It shows that when r is x dependent, S x does not 
follow exactly the spatial profile of S x . Due to this mis- 
match, according to Eq. (fl7|) , the S z component of the 
spin density is not zero on both sides of the transition 
layer and spin injection is feasible. 

The rest of the boundary conditions can be obtained 
from Eqs. (fl"5|) . (fT6|) . and (fl7|) by integrating these equa- 
tions across the layer and taking into account that S z 



and the sum of the first two terms in (fT8|) vary slowly 
within the layer. Therefore, the term proportional to 
\7 X S Z can be disregarded when integrating Eq. (fT5|) and 
also can be disregarded the sum of the second and third 
terms of the last line. They give rise only to corrections 
~ a 2 d/ A. Also, taking into account that outside the layer 
t 1 = a! = 0, we obtain 



[V X S* + 2m*aS z } X2 = [V X S X + 2m*aS z ] Xi , 
[V X S Z - 2m*a(S x - St)] X2 = [V x 5 2 - 2m*a(S x - S£)] X1 
V B S»| Xa = V x S y \ Xl . (21) 



Hence, we arrive at the boundary conditions correspond- 
ing to spin current conservation across the interface in its 
conventional definition (up to the order a 2 ), although in 
general such a current is not a conserved quantity. 
IV. CONCLUSIONS 



We have derived general diffusion equations for spin 
transport driven by electric fields in a quantum well sub- 
ject to Rashba spin-orbit interaction. The diffusion equa- 
tions are valid for spatially dependent mean free paths 
and spin-orbit interaction couplings. From these general 
diffusion equations, we can determine the boundary con- 
ditions for the spin-density and spin-current components 
in regions with different mean free paths and spin-orbit 
couplings via boundaries that are thicker than the mean 
free path, but thinner than the spin-precession length. 
We find that spin injection is possible when the electric 
field is perpendicular to the interface, and also when the 
electric field is parallel to the interface. In the latter case, 
spin injection is possible provided the boundary layer has 
a spatially dependent mobility. 
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